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Abstract 

The dynamics of globally coupled map lattices can be described in terms of 
a nonlinear Frobenius-Perron equation in the limit of large system size. This 
approach allows for an analytical computation of stationary states and their 
stability. The complete bifurcation behaviour of coupled tent maps near the 
chaotic band merging point is presented. Furthermore the time independent 
states of coupled logistic equations are analyzed. The bifurcation diagram of 
the uncoupled map carries over to the map lattice. The analytical results are 
supplemented with numerical simulations. 
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1 Introduction 



A lot of physical non equilibrium systems are governed by the motion of many relevant 
degrees of freedom. The most prominent example is hydrodynamic turbulence but 
also optical, magnetic, chemical and biological systems fall into this class (e.g [1] and 
references therein) . From the theoretical point of view it is often too difficult to solve 
the equations of motion directly even on the largest accessible computers. Hence one 
is forced to investigate mathematical model systems especially if one is interested in 
the pattern formation out of a random state. One class of such systems is given by 
coupled simple maps [2] . Several phenomena arising in the context of hydrodynamics, 
optics and solid state physics have been treated by such models [3, 4]. Although the 
actual relation between the basic equations of motion and the map lattice model 
is often not very well justified the simplicity of these models allow for a detailed 
numerical investigation. Furthermore rigorous approaches can be applied [5, 6, 7]. 
Hence it may be possible to understand the basic mechanisms of pattern formation 
from a random state. But the general theory is just at the beginning. For that reason 
it is desirable to investigate the dynamics of simple spatially extended model systems 
beyond numerical simulations. 

Frequently coupled map lattices with short range interaction are investigated. This 
kind of interaction models the diffusive character of the spatially extended system. 
But it seems to be too difficult to investigate theoretically. A much more simpler 
type is given by an all to all interaction. This coupling may be considered as a 
limiting case of a long range interaction [8] or as an approximation of a short ranged 
coupled system in a (spatially) high dimensional lattice [9]. Furthermore such type 
of coupling has been used to model e.g. the dynamics in multimode lasers [10], of 
Josephson junction arrays [11], or the dynamics in certain biological systems [12]. Its 
mean field like character greatly simplifies the theoretical approach [13, 14]. Hence I 
arrive at the model equation 

AT-l 

Here / denotes the single site map which should admit a chaotic motion, e the coupling 
constant, and N the systems size. The coupling is brought about by the mean field 

hn = f {xi^^) . (2) 

-'^ tj.=0 

Although much more general types of globally coupled systems may be considered I 
restrict the discussion to eq.(l). The limit of large system size iV 3> 1 is of special 
interest because in that case the number of degrees of freedom which are relevant for 
the time evolution may become large. Some surprising feature of global quantities like 
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the mean field (2) has been discovered recently [15]. The mean square deviation does 
not decrease with the system size but saturates at a finite value. This phenomenon 
has been termed "violation of the law of large numbers" and has been attributed to 
some hidden coherence in the coupled map lattice [16, 17, 18, 19, 20, 7]. 

The limit of large system size causes some problems. On one hand a direct numer- 
ical simulation of eq.(l) needs increasing computer capacities. On the other hand the 
occurrence of long transients considerably complicates the study of stationary states 
[21, 7]. For that reason an alternative approach is desirable. Because of the global 
coupling the systems allows for a description in terms of a reduced density [16, 22, 7]. 
Consider the probability that a phase space coordinate takes the value x. It reads 

1 Af-l 

pn{x) = ^ ^ - ^^:^) ■ (3) 

The mean field can be expressed in terms of this density via 

^n = j f{x)pn{x) dx . (4) 
For the time evolution of the density (3) one obtains from eq.(l) 

Pn+i(x) = y 5{x- Tn{y))pn{y) dy (5) 
where the mean field map T„ is given by 

Tn = F^of (6) 
= (l-e)x + e/i„ . (7) 

Eqs.(4), (5), (6) and (7) constitute a closed set of exact nonlinear evolution equations 
which is easier to handle as the original coupled map lattice. The system size enters in 
this formulation only via the structure of the density (3). In the limit of infinite system 
size the density tends for typical phase space points to a "sufficiently continuous" 
function^ (cf. ref. [23] for related rigorous statements). Therefore the limit of large 
system size is contained in this mean field description in a quite simple way. I restrict 
the analysis in the sequel to the discussion of eq.(5) and continuous densities. 

Let me focus on the investigation of stationary states, that means fixed points 
and periodic orbits of the mean field equation and their stability. For simplicity in 

^By this I mean that the corresponding measure Unig) '■= J g{x)pn{x) dx, V5, tends to a hmit 
which has a non atomic component. Although the formulation in terms of measures is physically 
less appealing one should keep in mind that equations like eq.(5) have to be understood in the weak 
sense, Hn+i{g) = MnCs °Tn), Vg, on a rigorous level. Especially the formal derivative p'n{x) has to 
be identified with the measure fnig) '■= — J g'{x)pn{x) dx, Vg, even if the derivative of the density 
is not defined. 
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the notation the exphcit formulation will be given for fixed points only. The general 
considerations are easily extended to periodic orbits by investigating the fixed point 
problem of a suitable iterate. The fixed points are determined by 



where denotes the mean field map evaluated at the fixed point p*. Without speci- 
fying the single site map less can be said about the structure of the solution. However 
some conclusions can be drawn about the stability. Consider a small deviation from 
the fixed point density pn{x) = p^{x) + 6pn{x). As long as linear stability analysis 
can be apphed these deviations obey (cf. eq.(5)) 



Stability is determined by the eigenvalue problem which corresponds to this linear 
evolution equation. Obviously this eigenvalue problem has the form of a perturbed 
Frobenius-Perron equation for the map T^. The second term in eq.(9) yields the 
formal perturbation. It is well known that the Frobenius-Perron equation has an 
eigenvalue one which is the largest in modulus. Because of a symmetry of the full 
evolution equation (5) and the normalization of the density this eigenvalue is persis- 
tent with respect to the perturbation [7]. The corresponding eigenfunction may be 
understood as a kind of Goldstone mode. If the map T^, has certain hyperbolic proper- 
ties then the largest eigenvalue is isolated. Hence one needs a finite coupling strength 
in order that additional eigenvalues cross the unit circle. It is therefore expected that 
stationary states are stable in hyperbolic systems at least for small coupling strength. 
On the contrary the spectrum of the Frobenius-Perron operator for non hyperbolic 
maps is usually degenerated on the unit circle [24, 25]. Therefore an infinitesimal 
coupling may induce instability and a complicated dynamical behaviour. 

On this level the considerations on stationary states and their stability is only 
qualitative. In section 2 the stability analysis will be put on a more rigorous basis. 
The explicit computation of stationary states requires the knowledge of the single 
site map /. In section 3 the analysis will be presented for the "hyperbolic" tent map 
which includes the bifurcations that occur in this simple system. The problem of 
non hyperbolic systems is analyzed in the context of the logistic equation. A partial 
bifurcation analysis and numerical simulations are presented in section 4. Finally I 
discuss the implications for the dynamics of the original coupled map lattice . 




(8) 




(10) 
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2 Stability of stationary solutions 



The stationary solutions depend on the special map lattice under consideration. But 
the stabihty can be discussed quite generally if the existence of the fixed point is 
presupposed. Therefore the computation of the stationary states is postponed to the 
subsequent sections. Here considerations concerning their stability are presented. 

I suppose that eq.(8) has been solved for a continuous density which should 
coincide with the Sinai-Ruelle-Bowen (SRB) measure of the map T*. The instability 
of such a distribution is indicated by an exponentially increasing solution of eqs.(9) 
and (10). Its formal solution reads^ 

5pn{x) = j 5{x-T:{y))5po{y)dy 

-Y^JT^.,.,^j6{x-T:{y)){T:)'{y)p4y)dy (11) 

57; = ej f{T:{y))p4y)dy 

""^ r d 

+6 STn-i-a / — / {T:{x)) p.{x) dx (12) 

(T=0 •' 

For any continuous initial condition 5pQ the first term of eqs.(ll) and (12) decays 
to zero because the SRB measure attracts any continuous distribution. Therefore 
an exponential increase, that means instability, can be excluded if for all sufficiently 
smooth functions g and h the quantity 

Jn ■■= J ^9 {T:{x)) ■ h{x)p4x) dx (13) 

remains bounded if n tends to infinity. The subsequent considerations focus on the 
discussion of this quantity. 

The expression (13) has a close relation to the linear response of the map 
[26, 27]. To clarify this point consider for the moment the map and denote its 
SRB density by p*. Let xq — x + eh{x) be an ensemble of initial points where the 
values X are distributed according to the SRB measure. The time evolution of the 
expectation value of a smooth observable g reads 

{g)n^{g{T:[x + eh{x)])). (14) 

where the brackets (. . .)* denote the SRB average. The quantity (14) tends towards 
the stationary value (51)* in the limit n ^ 00. Expansion in the small quantity e leads 
to 

{g)n-{g)*^e{^g{T:{x))-h{x)).-rO{e') . (15) 
^In the weak sense the relation (11) reads Sunig) = ^l^oig ° T^) — Sct=o ^T^n-i-al^*{{g ° T^)') . 
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The right hand side yields the formal linear relaxation function which coincides with 
the expression (13). As long as the expansion in e is uniformly valid in n this quantity 
decays because the left hand side does. But a linear response theory is in general not 
valid in the context of low dimensional maps even for hyperbolic systems. Never- 
theless any system which allows for the application of linear response theory leads 
to a bounded relaxation function (13) and therefore to stable stationary solutions in 
globally coupled maps. But this condition is too strong for general considerations. 

The estimation of the quantity (13) requires the knowledge of the density p*. It 
can be computed for the class of piecewise linear Markov maps. This kind of systems 
is characterized by maps of the interval / which are linear on a suitable partition 
h — [Oi^jQiz+i], < u < M, I — \J,^I„, — Tl(x) for x E Ii,. Furthermore it is 
required that the image of each interval is a union of other intervals T^{Iu) — U^/^. 
This property implies the expansiveness of the map. Although these systems seem 
to look rather special a lot of one dimensional maps can be described in terms of 
such models by considering fine partitions I^, of the phase space / [28, 25]. It is well 
known that the dynamics of such systems is equivalent to a subshift of finite type 
and that the invariant density is a piecewise constant function, p^{x) = p^, for x £ /j,. 
Integration by parts immediately yields for the relaxation function (13) 

„ M-l 

Jn^ 9 {T:{x)) h\x)p.{x) dx+Y.9 {T:{a,)) h (a,) (p, - p,_i) (16) 

The relaxation function is obviously a bounded function. Whereas the first contribu- 
tion decays to a constant the second term is periodic for large n. 

The considerations of this section imply that the stationary states of the mean 
field equation are stable if the mean field map T^, can be cast into the form of a 
piecewise linear Markov map. At this level nothing can be said about non hyperbolic 
systems like the frequently discussed logistic equation. In this case two different 
situations have to be considered. If the map T^, possesses a stable periodic orbit 
then the corresponding SRB density consists of a finite number of 6 peaks. But then 
the factors {T^)\x) in eqs.(ll) and (12) decay exponentially and imply the stability 
of the solution at least if the initial disturbance 5po is contained in the basin of 
attraction of the periodic orbit. This result is obviously trivial. In the other case 
the SRB density develops square root singularities'^. But these singularities induce 
an exponential increase of the relaxation function (13) [26] and turn the stationary 
solution unstable. This fact corresponds to the violation of linear response theory 
because of the structural unstable character of the map. Hence one might expect 
that a small coupling e induces a non stationary behaviour in the map lattice. We 
will come back to this phenomenon in section 4. 

^This statement is rigorous if the orbit of the critical point has finite length. 
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3 Stationary states for coupled tent maps 



A very simple but nontrivial model system is given by coupled tent maps 

f{x)^l-a\x\ . (17) 

The isolated map has a chaotic motion for 1 < a < 2 and shows a cascade of chaotic 
band merging at parameter values a„ = 2^/^". The mean field map (6) which governs 
the dynamics of the globally coupled map lattice reads 

Tn{x) = 1 — ea{\x\)n — (1 — e)a\x\ (18) 

where the brackets (. . .)„ denote the average with respect to the density Pn{x). 

For a dense set of parameter values a, e one can manage that the orbit of the 
extremal point a; = terminates at some unstable periodic orbit. The corresponding 
finite set of orbit points yields a Markov partition. The results of section 2 imply 
that any smooth periodic solution of the mean field equation is stable. Hence only 
the construction of these solutions from eq. (8) is desired. 

For that purpose it is useful to make a linear time dependent scale transformation 
X = '-fnZ in the mean field equation (5) in order that the mean field map (18) takes 
the form of the uncoupled map (17) . The evolution equation for the scaled density 

Pn{z) = hn\pn{lnZ) (19) 

reads 

Pn+l{z) - / ^ - Tn{z')) p{z') dz' (20) 

where the scaled mean field map is given by 

f„M = 1 - - N (21) 

if the scaling constants are determined by 

7„+i = 1 - ea|7„|(|2;|)^ . (22) 

In the sequel the brackets denote the average with respect to the scaled density (19). 
Fixed points: For time independent solutions p„ = eqs.(21) and (22) yield 

7* = l/(l + ea(|z|)*) > (23) 
f,{z) = l-(l-e)a|^| . (24) 

The actual shape of the fixed point density is determined by an equation analo- 
gous to eq.(8). It coincides with the smooth invariant density of the tent map at a 
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parameter value Qred. '■— (1 — The global coupling has the effect to reduce the 
parameter in the tent map. 

Period 2 solution: The single tent map (17) admits a period two chaotic solution 
for < a < \f2. A corresponding period two orbit p\iP\ may be expected in 
the mean field dynamics (20) too. It is determined by the following set of coupled 
equations 

~p\{z) = Js{z-f*{z'))pl{z')dz' 

= J S {z - (f* o f *)(/)) pliz') dz' (25) 

ft{z) = 1-^^N, 1-2 (26) 

72* = l-ea|7i*|(|^|)I, 1^2 (27) 

where (. . .)^/2 denote the average with respect to the densities 

First of all let me investigate the form of the general solution of eq.(25). For that 
purpose the ratio c := I2II1 considered as a free parameter which will be fixed 
later. The two times iterated mean field map which determines the shape of the 
density is given by 

{f*ofl){z) = l-ared.\c-ared.\z\\ . (28) 

Here 7* > has been assumed which is justified a posteriori. For a,.e,i. > \/2 the 
map admits only one continuous invariant density. It corresponds to the fixed point 
solution discovered above. Therefore let me concentrate on the opposite case. For 
< o,red. < and c e [1/ared., o^red] the map decomposes into two different ergodic 
components which are separated by the unstable fixed point (cf. Fig.l) <Fig.l 

The corresponding ergodic invariant densities are denoted by p'^^\z;c), \z\ < ({c) 
respectively p^'^\z; c), z > C{c). It is an elementary but important geometric property 
of the map (28) that the variation of the parameter c shifts the fixed point ({c) but 
leaves the two components of the map unchanged up to a linear scale transformation. 
Hence the dependence of the densities on the parameter c can be written as 

p'-\z;c) = ^p(-)f^.;c'V |.|<C(c) (30) 

-^);c'V z>ac) . (31) 
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In addition both densities are not independent of each other. The geometric properties 
of the map (26) imply that (appendix A) 

Js(x- f*{y)) p^^\y- c) dy = p(^) (x; (32) 

holds. The general solution of eq.(25) is an arbitrary convex combination of both 
ergodic components 

pI{z) ^ {l-a)p^-\z;c)+ap^+\z;c), a e [0, 1] . (33) 

But the parameter c is not at our disposal. It has to satisfy the additional con- 
straints (27). Because one has control of the c dependence of the solution via eqs.(30) 
and (31) it will be shown that the relative weight a of the two ergodic components 
determines this ratio. To this end let me first mention that thanks to eqs.(32) and 
(33) the relations 

(l^l)t = (l-«)(k|)(-)(c) + «(|.|)W(c) (34) 
{\z\)* = (l-a)(|;.|)W(l/c)+a(|;.|)(-)(l/c) (35) 

hold where the abbreviations 

(|^|)W(c):= J\z\p^^\z;c)dz (36) 

have been used. Then eqs.(27) read 

1/71* = c + ea[{l-a){\z\y-\c)+a{\z\y^\c)] (37) 
1/7; = l/c + 6a[(l-a)(|z|)W(l/c) + a(|^|)(-)(l/c)] (38) 

which result in 

l-c = ea[{l-a) {{\z\y-\c) - c{\z\)^^\l/c)) + a ((k|)W(c) - c{\z\)^-\l/c) 



(39) 

Because of the scahng property (30) and (31) the dependence of the expectation 
values on the parameter c can be evaluated to 

m^^\c) = l-i^(l-(|.|)(+)(l)) . (41) 



i-c(i: 

Here c' has been set to unity for convenience. If one additionally uses the relation 
(|z|)*-"'"''(l) = 1 — ared.{\z\)^~\l) between both expectation values (cf. appendix A) 
one ends up with 

^^l^ ea(ared. - l)(l-2a)r/ 

ared. - 1 + ea [1 - ttred.V + avictred. " 1)] 

9 



where the abbreviation rj — 1 — [1 + ared.]{\z\)''~\l) has been introduced. The right 
hand side depends on « in a monotonous way (cf. appendix A). Therefore for a 
values in a neighbourhood of a = 1/2 eq.(42) determines the unique c value in the 
allowed interval [l/ared.,cired.] which satisfies eq.(27) by construction. Hence every 
density of the form (33) yields a period two orbit of the mean field equation (20). 
The fixed point solution which also exists at the parameter values under consideration 
is contained in this continuous family for a particular value of the parameter a 

Higher order periodic orbits: It is obvious that the computations made above can 
be applied to any band merging in the tent map (e.g. consider the ergodic component 
p'--^^). Therefore at the critical parameter values (1 — e)a = 2^/^" a continuum of 
period 2" orbits occurs. Every solution is stable. Depending on the initial condition 
a particular state is attained in the course of the dynamics. 



4 Analysis of coupled logistic equations 

The most popular system which is discussed in the literature is given by coupled 
logistic equations 

f{x) = 1 - ax^ . (43) 

Although the dynamics of the uncoupled map is fairly well understood on a rigorous 
level [29] it is surprisingly difficult to obtain related results for the mean field coupled 
map lattice. The main problem originates from the non hyperbolicity of the map 
which causes the dense set of windows of stable periodic orbits in the uncoupled 
case. Because of the lack of suitable analytical approaches one is forced to resort to 
numerical solutions of eq.(5). Nevertheless some statements on stationary states can 
be made. 

Stationary states: The mean field map which governs the dynamics of eq.(5) reads 

T^{x) = 1 - ea{x^)n - (1 - e)ax'^ . (44) 

To analyze the stationary states it is again useful to introduce rescaled phase space 
variables so that eq.(44) takes the form of the logistic equation. Choosing 

7n+i = 1 - eajl{z'^)n (45) 
the mean field equation is cast into the form (19), (20) with 

f{z) = 1 - "--^7; (46) 

where the brackets (. . .)„ denote the average with respect to the density (19). Let me 
focus on the discussion of fixed points. Periodic orbits can be treated in principle by 



10 



considering the iterated system. However the discussion of the fixed point problem 
seems to contain the main structure of the whole bifurcation diagram. The fixed point 
density is determined by an equation analogous to eq.(8) where 

f^z) ^ 1 - ared.l*z^ (47) 

denotes the mean field map and the scaling parameter obeys 

7. = l-ea7,'(z'). . (48) 

If one regards for the moment the quantity c := ared.1* ^ ^^e parameter then 
is the unique SRB measure of an ordinary logistic equation^. This density fixes the 
quantity (z^)* — {z^)*{c). If one now eliminates the scaling parameter 7* from eq.(48) 
in favour of the abbreviation c one ends up with 

{l-eya = c + ec{c{z^),{c)-l), c G [0, 2] . (49) 

Eq.(49) yields a family of curves in the parameter plane (a, e). On every curve the 
coupled map lattice admits a stationary solution. If the parameter c is chosen in such 
a way that its value corresponds to a periodic window in the logistic equation then 
the density consists of a finite collection of 6 peaks. Due to to the discussion at the 
end of section 2 these states are dynamically stable. In the opposite case less can be 
said about the density p^, and their stability. However the remarks made in section 
2 indicate that the solutions are dynamically unstable at least for situations which 
correspond to Smale complete logistic equations. 

To sum it up it can be said that eq.(49) provides a foliation which translates 
the bifurcation structure of the logistic map to the parameter plane (a,e). States 
which correspond to periodic windows lead to dynamical stable fixed points. The 
fixed point analysis presented here has much in common with the approach used in 
ref. [18]. But in contrast to this approach eq.(49) provides an explicit expression for 
the bifurcation lines and clarifies the stability of the solutions. Finally it should be 
mentioned that to my best knowledge it is not clear whether the parametrization 
(49) depends continuously on c^. Additionally it is not clear whether the lines can 
intersect so that coexisting stationary states occur. 

Numerical simulations: The discussion in the in the preceding paragraph has 
shown that the bifurcation diagram becomes tremendous complicated even if only 
the fixed points are considered. Although there exist stable solutions it is not clear 

^Although every physicist believes that the logistic equation possesses an SRB measure I am not 
aware of any proof that eq.(43) possesses for every a G [0,2] a unique measure that describes the 
evolution of Lebesgue almost all initial points x. 

^Whether (2;^)*(c) depends continuously on c seems to be an unsolved problem. In fact both 
hypothesis can be found in the literature[18, 22]. 
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whether a typical initial condition settles on these solutions. But it is expected that an 
infinitesimal coupling changes the dynamical behaviour in contrast to the hyperbolic 
example treated in the preceding sections. 

To shed some light on these problems the evolution equation (5) has been iterated 
numerically. The numerical algorithm which has been used is briefly described in 
appendix B. Main interest has been focused on the regime of weak and moderate 
coupling strength. The uniform distribution was chosen as a typical initial density. 
It corresponds to a random uniform distribution of phase space coordinates for the 
original coupled map lattice (1). The results described below seem not to depend 
on this choice. If the parameter a is restricted to a large periodic window one finds 
in the weak coupling regime stable periodic solutions of the same period. They 
correspond to the solutions which have been predicted by the analytical approach 
presented above. More interesting things happen if the parameter is chosen as a 
"typical" chaotic value. Several results for the choice a = 2 are presented here but 
qualitative similar results are obtained at all other values that have been investigated. 
Even for small coupling strength the system does not settle to a stationary state. The 
distribution of the mean field develops a complicated shape which shows an apparent 
nonsystematic variation with the coupling strength (cf. Fig. 2). The width of the <Fig.2 
distribution is finite and is responsible for the so called " violation of the law of large 
numbers" [22, 7]. The time evolution of the mean field is conveniently analyzed in 
terms of its power spectrum (cf. Fig. 3). For weak coupling strength the spectrum <Fig.3 
shows a small noise level on which sharp peaks are superimposed. By increasing the 
coupling strength the peaks broaden and the noise level increases. The actual shape 
of the spectrum seems to depend sensitively on the coupling strength. Because of the 
comphcated structure of the partial bifurcation diagram described above it seems to 
be hopeless to attribute the peaks to definite periodic orbits. Finally Fig.4 contains <Fig.4 
the time cvohition of the density p„ for two values of the coupling strength. In the 
case of weak coupling the density looks like a stochastic perturbation of the invariant 
density of the uncoupled system. The case of moderate coupling for which a broad 
band noise in the power spectrum of the mean field is observed yields a density with 
a strongly developed structure. It possesses sharp peaks which show an intermittent 
time evolution. The appearance of localized peaks in the density indicate a clustered 
state in the original coupled map lattice. Hence this type of motion can be viewed as 
an intermittent dynamics of clusters in the original system. Qualitatively this type of 
motion can be attributed to the stable period three window in the logistic equation. 
Indeed for slightly different coupling strength a stable period three state can be found 
in numerical simulations of the mean field equation. 
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5 Conclusion 



The simplicity of the global coupling has allowed for a mean field like description of 
the coupled map lattice. It contains the limit of infinite system size in a quite simple 
way namely by a "smoothing" process of the reduced density. In a certain sense this 
description circumvents the problem of supertransients. If the transient time increases 
tremendous with the system size, e.g. exponentially, then the relevant dynamics of 
the spatially extended system may be the transient one and not the (mathematical) 
stationary state. This behaviour is of course difficult to analyze analytically for the 
original system (1). One merit of the mean field formulation consists in the fact that 
in the limit of infinite system size, that means among continuous densities, these 
transients correspond to stationary states which can be discussed more easily. 

For the example of tent maps a fairly complete survey over the stationary states 
has been given. It is worth to mention that these solutions can be found in numerical 
solutions of the mean field equation as well as direct simulations of the coupled map 
lattice. They are attained for all initial conditions so that these states are appar- 
ently not only locally but also globally stable. Hence the stationary state is a unique 
fixed point for 2 > (1 — e)a > \pl. The finite size effects induce fluctuations around 
this flxed point so that global quantities like the mean fleld "obey the law of large 
numbers". At other parameter values one solution is selected from the continuum 
of periodic states depending on the chosen initial condition. Because of its period- 
icity the mean square deviation of global quantities saturates at a flnitc value and 
apparently "violates the law of large numbers". The value of saturation depends on 
the initial condition (cf. [16]). At different parameter values that are not covered 
by the present approach, especially at negative coupling, additional bifurcations may 
occur [22]. These values have been skipped in the present discussion because the map 
lattice may have diverging trajectories. 

Beyond the phenomena encountered in the hyperbolic tent map the non hyper- 
bolic logistic equation shows additional features. First of all the degeneracy of the 
spectrum of the transfer operator causes, that even an infinitesimal coupling changes 
the dynamics drastically. This is obvious from the bifurcation analysis of fixed points. 
In contrast to the hyperbolic situations the numerical simulations indicate that the 
locally stable solutions are not attained for a typical initial condition. Hence their 
domain of attraction seems to be small. In view of the complicated structure of 
the bifurcation diagram with probably infinitely many coexisting stable states this 
observation is not very astonishing. The lack of a unique stable fixed point also ex- 
plains that the phenomenon of the "violation of the law of large numbers" is typically 
observed in non hyperbolic examples at small coupling strength. 

Finally it should be mentioned that it is difficult to decide which features survive 
if the condition of infinite coupling range is relaxed. It is however remarkable that 
hyperbolic coupled maps show often a rather simple evolution whereas non hyperbolic 
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maps show a complicated dynamics with many coexisting stable or metastable states 
and intermittent behaviour. The investigation of globally coupled systems may shed 
additional light on these phenomena. But further investigations are required. 
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Appendix A 

Proof of eq. (32): The abbreviation 

^(^;l/c):=l-^k| (50) 
is introduced to take explicitly the c dependence into account. Obviously the relations 

f*{z)^g{z-l/c), f*{z)^g{z-c) (51) 

hold. The ergodic components p^^\z] c) are the continuous invariant measures of the 
map g{g{z]l/c)]c). Hence the functions 

^P^^\z) := j 6{z-g{z'-Mc))p^^\z';c)dz' (52) 

are continuous invariant densities of the map g{g{z; c); 1/c) which coincides with the 
former one if c is replaced by 1/c. This map has two ergodic components which are 
contained in the two intervals [— C(l/c), C(l/c)] respectively [(^(l/c),l]. The original 
densities p^^^ are contained in the intervals [— C(c),C(c)] respectively [C(c),l]. One 
easily calculates that the images of these intervals obey 

^([-C(c),C(c)];l/c) = [C(l/c),l], ^([C(c),l];l/c) C [-C(l/c),C(l/c)] (53) 

Hence the densities (52) are continuous functions on the intervals (53) and coincide 

for that reason with the ergodic components p^^\z; 1/c). 

Relation between {\z\)^^^ and {\z\)^'^^: The case c = 1 will be considered throughout 
this paragraph and the argument c = 1 will be suppressed in the notation. On the 
two different ergodic components the map (28) which is simply the second iterate of 
the tent map reads 

g^~\z) := 1- aredX^- ared.\z\), \z\<C (54) 
g'^^\z) := 1 - ared.l'i- - dredM^ C < -^^ < 1 ■ (55) 
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Both maps are conjugate to each other via g'^'^^o/i = /iogr( ^ where h{z) — {l — z)/ared.- 
On one hand the conjugacy imphes 

J zp^-^ {z) dz = J h-\z)p^^^ {z) dz=l- ttred. J zp^^^ {z) dz . (56) 
On the other hand the invariance of the density yields 

I zp^-\z) dz^l g^-\z)p^-\z) dz^l- area. + a^.d^l)^"^ ■ (57) 
If one combines both equations one obtains 

(|z|)(+) = l-a,erf.(k|)(-) . (58) 

Discussion of eq.(42): It will be shown that eq.(42) determines a unique non 
negative value for every admissible a value. Because p^'^z;!) is contained in the 
interval [-((!),< C(l)] the inequality (|2|)(-)(1) < C(l) = l/ia-red. + 1) holds. But 
then < 77 < 1 is valid. The denominator of eq.(42) obeys for a e [0, 1] 

dred. - 1 + ea{l - ttred.) < dred. - 1 + ea[l - ared.V + C^Vi^^red. " 1)] ■ (59) 

Hence it is positive for e < 1/2 (this bound on the coupling strength can be relaxed if 
a better estimate of the quantity t] would be derived) . But then the right hand side 
of eq.(42) is a monotonous function so that for every ck e [0, 1] a unique value of c is 
found. This value is positive because 

^ ^ _ ea{ared. - l)v 



O-red. - 1 + ea[l - ared.1l + ^(Ored. " 1)] 

{ttred. - 1)(1 - ea) 



ared. - 1 + ett[l - ttred.V + Vi^ired. " 1)] 

holds for e < 1/2. From the estimate 

^a{ared. - 1)(1 - 2a)?7 



> (60) 



dred. - 1 + ea[l - ttred.V + Oir){ttred. - 1)] 



< (61) 



ea 



one additionally obtains that for sufficiently small coupling strength e any value a e 
[0, 1] corresponds to a unique value c e [1/ ttred., O'red.]- 



Appendix B 

Let me briefly describe the basic idea for the numerical simulation of eq.(5). The 
straightforward approach consists in an equidistant partition of the interval and the 
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approximation of the density by a piecewise constant function in every iteration step. 
But an equidistant partition seems to be not very suitable especially if the density 
develops square root like singularities. A finer partition in those regions is desirable 
whereas in flat regions a coarser partition is sufficient. For that reason the densities 
are approximated by piecewise constant functions so that every interval contains the 
same weight. To be definite suppose that at time n a partition 7^"^ = [^^1-^)^1"^]) 
1 < k < N in N intervals is given so that the density has the form 

N ^ 

Pni^) = E (n) A:jW (^) (62) 

k=i NAal ^ 

where Aa^""* := a^""* — a^k-i denote the lengths of the intervals and xj stands for the 
characteristic function of the interval J. The partition at time n + 1 is determined in 
such a way that the density yields the same weight for every interval. This 
prescription results in 



^ = [i,, lsiz-Uz'))p^{z')dz'dz 

= y [,\ // 6 (z - TJz')) dz' dz ^ (63) 



The partition at time n + 1, {a^k'^^^^}, can be determined easily from eq.(63) in an 
iterative way because the remaining integrals are nothing else but the lengths of the 
intervals Ii^^ fl T~^(ljr'^^^)^. At time n + 1 the density is again approximated by 
a piecewise constant function on the new partition (cf. eq.(62)) which corresponds 
formally to the canonical orthogonal projection. 

Numerical estimates show that the error of the algorithm is of the order 0{1/N) 
even if the density develops singularities. For the numerical simulations a partition 
of size N — 2000 has turned out to be sufficient so that the computations can be 
performed on every PC. 



^An explicit algorithm for symmetric unimodal maps is available on request. 
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Figure captions 



Fig.l Diagrammatic view of the two times iterated mean field map (28) for \fl < 
dred. < \/2- The diagonal, the dashed boxes, and the unstable fixed point C(c) 
are indicated for clearness. 

Fig. 2 Distribution function of {x^)n = (1 ~ hn)/a for coupled logistic equations at 
a — 2 and several coupling strengths. The distributions have been obtained 
from a time series of 10^ data points. They are displayed as histograms with 
200 boxes on the abscissa. The size of the boxes are adapted to the width of 
the distribution. 

Fig.3 Power spectra of for coupled logistic equations at a = 2 and several 

coupling strengths .The spectra have been obtained from the Fourier transform 
of a time series of length 1024. An average over 200 different series has been 
performed. 

Fig.4 Time evolution of the reduced density Pn{x) for coupled logistic equations at 
a = 2 and two values of the coupling strength (a) e = 0.01, (b) e = 0.2. The 
densities are displayed as histograms with 100 boxes on the x axis. 
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